Vibrating soap films: An analog for quantum chaos on billiards 
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We present an experimental setup based on the normal modes of vibrating soap films which shows 
quantum features of integrable and chaotic billiards. In particular, we obtain the so-called scars 
-narrow linear regions with high probability along classical periodic orbits- for the classically chaotic 
billiards. We show that these scars are also visible at low frequencies. Finally, we suggest some 
applications of our experimental setup in other related two-dimensional wave phenomena. 



I. INTRODUCTION 

In recent years, there has been increasing interest 
in the properties of quantum systems whose classical 
analogs are chaotic. ||f|] Part of the work in this new 
field, called quantum chaos, [Q refers to essentially two- 
dimensional cavities or wells of infinite potential called 
billiards. These billiards can take different forms, such as 
rectangles, circles and other more complicated geometries 
(see Fig. |l|). The circle and the rectangle correspond to 
integrable systems. Q Furthermore, we include in Fig. |l| 
the so-called Bunimovich stadium and the Sinai billiard 
which are completely chaotic. Q However, intermediate 
situations, i.e., systems with both integrable and chaotic 
behaviors, are the most common type of dynamical sys- 
tems, u 

The quantum analog of a classical billiard is called 
a quantum billiard and its eigenfunctions are closely re- 
lated to the classical features of the billiard. (?]] Quantum 
billiards obey the Hclmholtz equation with vanishing am- 
plitude on the border (homogeneous Dirichlet boundary 
conditions) . Such systems can be simulated as a drum or 
any other membrane vibrating in a frame. The principal 
purpose of this paper is to show that the quantum behav- 
ior of classically chaotic and integrable billiards can be 
modeled in a classroom with an analog experiment. The 
experimental setup basically contains a function genera- 
tor and a mechanical vibrator and is based on the vibra- 
tions of a soap film. || Thus, as Feynman said in 1963, 
"the same equations have the same solutions" . || 

In the next section we briefly discuss classical and 
quantum billiards. In Sec. Ill we introduce our exper- 
imental setup and show the analogy with the quantum 
billiard. Other uses of our experimental setup are dis- 
cussed in the same section. Some remarks are given in 
the conclusion. 

II. CLASSICAL AND QUANTUM BILLIARDS 

In order to make a more explicit definition of the clas- 
sical billiard we take a two-dimensional region denoted 
by R and define the potential V for the particle as 



V = 







I? 



oo otherwise. 



(1) 



boundary, it bounces following the law of reflection. Un- 
der this dynamics, the rectangle and circle billiards are 
regular. Typical trajectories inside are shown in Fig. |^. 
On the other hand, the dynamics of a particle in the 
stadium as well as in the Sinai billiard, are chaotic. Al- 
most all trajectories for these billiards are ergodic and 
exponentially divergent. Roughly speaking, this is so be- 
cause in the Sinai billiard, two very close trajectories are 
separated when one collides (or both) with the central 
circle. After some time, the separation between the tra- 
jectories (in phase space) is exponential. In the stadium 
two particles with very close initial conditions are focused 
when they collide with one semicircle. After this focus- 
ing, they began to separate until they bounce again but 
now on the other semicircle. The exponential divergence 
appears because the separation time is greater than the 
focusing time. Apart from these trajectories there also 
exist periodic orbits. These trajectories are unstable and 
typically isolated. Their number increases exponentially 
as function of their length, but they are of measure zero 
in phase space. We may find also families of unstable 
and non-isolated periodic orbits such as the "bouncing- 
ball" orbits in which the particle bounces between the 
two parallel segments of stadium. Some periodic orbits 
for the Bunimovich stadium are shown in Fig. |^. 

The time independent Schrodinger equation for the 
potential defined in Eq. (1) is 



V 2 * 



fc 2 * = in R, 

$ = on the boundary of R, 



(2) 



This means that inside R the particle is free and moves 
in straight lines. When the particle collides with the 



with the wave number k = (2mE/h) 1 ^ 2 . Here, E and 
m are the energy and mass of the particle and h is the 
Planck constant. The homogeneous Dirichlet boundary 
condition is obtained because if V = oo the wave function 
vanishes. The Hclmholtz equation in R and the Dirichlet 
boundary condition define the quantum billiard. Note 
that this is just the equation for the normal modes of 
a membrane if we interpret the functions as vibration 
amplitudes. 

For classically integrable billiards, the eigenfunctions 
are well-known. For example the solutions for circular 
and rectangular billiards are Bessel functions and sinu- 
soidal functions, respectively. |nj On the other hand, 
the features of wave functions for classically chaotic bil- 
liards have been well studied numcr ically by Heller. § 
Recently, experiments in microwave cavities have been 
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performed, jllj In Fig. |i| we show eigenfunctions of the 
Bvmimovich stadium we calculated numerically using the 
finite element method. |12| However they can alterna- 
tively be calculated by standard software. The eigenfunc- 
tions of figures_0(a)-(c) show certain similarities with the 
orbits of Fig. yTFoHowing Heller, we say that the eigen- 
functions are "scarred" by the orbits. Fig. ||(d) shows 
what is called a "whispering gallery" state, because there 
exist certain galleries [T^ ] in which the sound travels in- 
side them, following the border. This kind of state is 
associated with orbits also close to the boundary. The 
eigenfunction shown on Fig. ||(e) resembles noise Q when 
we see only a quarter of stadium. 

The appearance of the scars is quite well understood 
Jl4j based on the theoretical work by Selberg, Gutzwiller 
and Balian. jl5| Due to the low density of the short pe- 
riodic orbits they may well be seen in quantum experi- 
ments and simulations cither as dominant features in the 
Fourier spectrum or as scars. While the Fourier analysis 
of experimental data in atomic |l6| and molecular physics 
jl7j is quite striking, direct demonstrations of scars are 
difficult even in microwave cavities. A simple explana- 
tion of scarring is based on de Broglie waves. Close to 
the periodic orbit there exist standing de Broglie waves 
whose wavelength A is associated to the length L of the 
periodic orbit: 



2L = n\ 1 7i=l,2,3, . 



(3) 



These de Broglie waves are localized around periodic or- 
bits due to the exponential divergence of nearing tra- 
jectories. In the next section we will show how a simple 
demonstration setup can display the most interesting fea- 
tures of the eigenfunctions on a soap film. 

III. SOAP FILM ANALOGY 



The normal modes of a rectangular and circular soap 
film have been well studied. []18|— p0[| A textbook showing 
these eigenfunctions is French's book entitled Vibrations 
and Waves, jllj The governing equation is the time inde- 
pendent wave equation (Eq.||), but now with a different 
interpretation: '5 is the membrane vibration amplitude 
and the wave number is now k — lj/v, with lu the angular 
frequency and v the speed of the transverse waves on the 
membrane. 

A problem arises for the experimental setup of this 
analog: At high frequencies (corresponding to the semi- 
classical limit), the damping is large. To solve this prob- 
lem we feed energy into the system permanently with 
an external resonator of a well-defined but variable fre- 
quency. We chose to feed the external frequency into the 
system by vibrating the wire that delimits our soap film. 
We use a mechanical vibrator (PASCO Scientific model 
SF-9324, see Fig. ||) connected to a function generator 
(we used a generator Wavetek model 180 [^2[) and used 
wires with different shapes (rectangle, circle, Bunimovich 
stadium, Sinai billiard, and some other of interest |23|) 
Alternatively a speaker could be used to transmit the fre- 
quency through the air. |l^,^,^| The chemical formula 
for a soap film with large duration is given by Walker, 



p3 but can be made up with soap, water and glycerin 
by trial and error. 

We can now start the demonstration. In Fig. |](a) we 
show a normal mode of the rectangle, in agreement with 
the known result. The Fig. ^(b) shows a normal mode for 
a circle displaying a Bessel function. Roughly speaking, 
the shining and dark zones establish a periodic pattern 
associated to the normal mode. Although the pattern es- 
tablished cannot give a quantitative measure of the am- 
plitude, it is sufficient to give an idea of the form of the 
normal mode and to identify it. However, the more inter- 
esting normal modes are some of the classically chaotic 
billiards. In Fig. |6|(c) we show a "bouncing-ball" state at 
low energies for a Bunimovich stadium. The alternating 
dark and shining zones in this and in the following figures 
make evident the presence of standing waves in the mem- 
brane. These waves are associated to de Broglie waves in 
the quantum billiard and at the same time they are as- 
sociated with periodic orbits in the classical billiard. In 
Figs. ||(d) and (e) scarred eigenfunctions are displayed, 
the latter in the high-frequency regime. In order to show 
that they are scars and not some spurious effect of our 
very simple experimental setup, we calculate numerically 
the normal modes at frequencies near the experimental 
ones. We observe in Fig. | some of these eigenfunctions 
and the corresponding orbits in Fig. ||. 

We want to mention that the normal modes in soap 
films can also be used for other two-dimensional wave 
phenomena, such as the search of normal modes of the 
clay layer corresponding to the old Tenochtitlan lake, 
which plays a crucial role in the earthquake damage pat- 
terns of Mexico City. ]2j| The application in this case 
comes from the fact that the upper clay layer of the 
Mexico Valley, as well as the soap films, are practically 
two-dimensional. J2q| A wire may readily be shaped to 
the corresponding boundary and the results are shown in 
Fig. |(f). 

Moreover, the experimental setup which we present 
may be used to show other related wave phenomena. 
Typical examples are scars on liquids, |Q Faraday waves 
-crystallographic patterns in large- amplitude waves ]2^]~ 
, cuasi-crystalline patterns on liquids, |2q ] sand dynamics 
p9[ or normal modes of Chladni's plates. |3^] These can 
be done by replacing the wire in our experimental setup 
with water containers or by thin plates. If we want to see 
scars on surface waves, we must put a stadium-shaped 
tank filled by water. Faraday waves can be obtained 
adding shampoo to the water and increasing the ampli- 
tude and frequency but decreasing the level of the liquid 
up to several millimeters. It is not necessary to use tanks 
with different shapes because the patterns do not depend 
on the boundary. If we want to observe quasi-crystalline 
patterns with this experimental setup, we must excite the 
mechanical vibrator with- at least -two frequencies. pSj ] 
This is easily obtained by changing the sinusoidal time 
dependence of the of the driving force to a triangular 
one. |3l]] Another application is the demonstration of the 
modes of thin plates. In this case, as well as in the wa- 
ter tanks, the whispering gallery states are easily visible. 
We may use ellipses, pentagons or any other shape. All 
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these plates are excited in a point in which the mechani- 
cal vibrator loads them. Figure 0(a) shows a whispering 
gallery state for a stadium-shaped iron plate. Fig. ^(b) 
shows a scarred pattern for the same plate. Finally, if we 
consider a tank with a layer of sand of variable thickness, 
we may study several topics on dynamics of granular me- 
dia. As an example we can study "standing waves" on 
sand. jH 

IV. CONCLUSIONS 

The analog model of soap films for the quantum bil- 
liard gives a demonstration of quantum chaos features. 
For integrable regions the normal modes correspond to 
the eigenfunctions of two-dimensional square box and cir- 
cle. For classically chaotic billiards the normal modes 
correspond to scarred eigenfunctions in the scmiclassical 
limit. The experimental setup presented here is particu- 
larly cheap, simple and elegant- i.e. the alternative expe- 
riments do not satisfy the Helmholtz equation and/or 
boundary condition. Microwave cavities are expensive 
and do not display the scars in a directly visible ma- 
nner-. Thus the simplicity of the experimental setup 
and the facility to put any shape, makes it very suit- 
able for the undergraduate laboratory. Furthermore the 
normal modes of soap films can be used to demonstrate 
other two-dimensional analog phenomena. Finally, we 
want to mention that the proposed experimental setup 
can be quickly changed to study other related wave phe- 
nomena. A wide variety of highly nontrivial wave-like 
phenomena can be displayed with minimal experimental 
requirements. 
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FIG. 1. Integrable billiards: (a) rectangular and (b) circu- 
lar. Chaotic billiards: (c) the Bunimovich stadium, formed 
by two semicircumferences of radius r joined by two segments 
of length d; (d) the Sinai' billiard, formed by a square with a 
circle of radius r inside 




FIG. 2. Typical trajectories in billiards: (a) in the rectan- 
gle, (b) in the circle (This trajectory shows a caustic), and (c) 
in the Bunimovich stadium. 
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FIG. 3. Periodic orbits in the Bunimovich stadium. Notice 
that in some cases more than one orbit is present and that all 
obey the reflection law. 
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FIG. 4. Eigenfunctions of the Bunimovich stadium calcu- 
lated by the finite element method for a quadrant of stadium 
with Dirichlet boundary conditions. The maxima and minima 
of each normal mode correspond to the darkest zones; (a) a 
typical state scarred by a short classical orbit connecting the 
two extremes of the stadium, (b) a state scarred by an orbit of 
large period, (c) a typical state scarred by a "bouncing-ball" 
orbit, (d) "whispering gallery" state, and (e) a typical noisy 
state if one looks only at a quadrant. 



FIG. 6. Normal modes for different shapes, (a) Low mode 
for the rectangle at a frequency of 22.48 Hz (The frequency 
changes slightly for different soap solutions and different 
thicknesses of the soap film.) The rectangle is of size 0.215 
mx 0.155 m. (b) Bessel function for a circular-shape wire. 
The radius of the wire is 0.092 and the excitation frequency 
is around 33.5 Hz. (c) A "bouncing-ball state" near to 40 
Hz. (d) A scarred state corresponding to the numerical one 
showed in Fig. ^(b) at frequency close to 110 Hz. (e) A scarred 
state by a short classical orbit connecting the two extremes 
of the stadium for a Bunimovich stadium with r — 0.035 and 
d = 0.035 . The corresponding frequency is about 478 Hz. (f) 
A normal mode for a wire with "old Tenochtitlan lake" form 
at frequency v = 48.88 Hz. 




FIG. 5. Experimental setup used to show the quantum 
features of chaos. Notice that the membrane is excited in a 
way that breaks the symmetry. 




FIG. 7. Nodal patterns observed on iron plates with sta- 
dium shape. The plate was of 1 mm thickness and we used 
r = 0.1m and d = 0.1 m. (a)A whispering gallery at 5.79kHz. 
Notice that this pattern is better defined near the boundary 
and between the two segments of the stadium. This means 
that the amplitudes are higher in these regions, (b) A scar 
which reflects the orbit connecting the two extremes of the 
stadium. In this case the pattern is better defined on the 
periodic orbit. The frequency for this case was around 4.22 
kHz. 
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